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[Text] 


Abstract 


An experimental impulse radar based on the Transient 
Electromagnetic Measurement System (TEMMS) has 
been used to measure the transient electromagnetic (EM) 
pulse responses of a conducting sphere and a conducting 
prolate spheroid in actual space; this technique is not 
subject to the limitation of the traditional Time-Domain 
Scattering-Range Technology. The impulse responses of 
the conducting sphere and the conducting prolate 
spheroid are derived from the measured transient pulse 
responses using the modified steepest descent method. A 
transient response method for measuring the wideband 
radar cross section (RCS) of targets in the time domain is 
presented. This method is used to measure the RCS and 
resonance frequencies of the conducting sphere and the 
conducting prolate spheroid. 


I. Introduction 


The transient EM pulse signal has an extremely wide 
bandwidth which extends from very low frequency up to 
microwave frequency; it is an electromagnetic source 
that has not been effectively exploited by mankind. For 
example, in target detection and identification, the tran- 
sient EM pulse signal can be used to extract a wide 
variety of target information. It is a signal that has great 
potential value in high-accuracy measurement, classifi- 
cation, identification and imaging of targets in free space 
(e.g., airplanes, missiles, space vehicles and decoys) as 
well as targets in dissipative medium (e.g., underground 
cables, cavities and land mines). The transient target EM 
characteristics also provide the basis for ultra-wideband 
(UWB) impulse radar counter-stealth technology. 


To study the transient target EM characteristics, it is first 
necessary to obtain the transient EM response; there are 
two approaches for obtaining transient EM response: the 
frequency-domain approach and the time-domain 
approach. 


In the frequency-domain approach, a single-tone signal 
is used to illuminate the target, and the amplitude and 
phase of the return signal are measured. In order to 


eliminate the effects of surrounding disturbances, the 
measurements must be performed in a microwave ane- 
choic chamber. Furthermore, in order to obtain the 
transient response, measurements must be made at suc- 
cessive points along the entire frequency band that 
contains the frequency spectrum of the signal. For 
example, at the Electronic Science Laboratory (ESL) of 
Ohio State University,' 1,199 amplitude and phase 
measurements of the return from a 6-inch conducting 
sphere are taken at 10-MHz increments in order to 
obtain its impulse response; the amplitude and phase 
data are processed using the discrete inverse Fourier 
Transform to obtain the time-domain waveform. 
Because the spectrum is abruptly cut off at 11.99 GHz, 
the waveform contains obvious advance waves (i.e., 
waves that appear before the signal reaches the target) 
and high-frequency oscillations that have no physical 
meaning; they must be removed by applying a window 
function to obtain an approximate impulse response. 
Therefore, the frequency-domain approach involves 
many complications. 


The conventional time-domain approach uses the Time- 
Domain Scattering-Range Technology (TDSRT),*~* 
where a large conducting plate is set up horizontally in an 
indoor facility (Figure 1); the measuring instrument and 
equipment are placed underneath the conducting plate. 
and the monocone (monopole) transmitting antenna is 
placed vertically on the plate. One-half of a mirror- 
symmetric target (e.g., a conducting sphere) is placed on 
top of the plate, and a monopole receiving antenna is 
placed close to the target. The pulse signal 1s transmitted 
by the monocone antenna; the direct signal and the 
scattered signal from the target are received by the 
monopole antenna and sent to the sampling oscilloscope. 


TDSRT is based on the mirror-image principle; it is 
simple and requires no microwave anechoic chamber. It 
can measure the transient EM pulse response of mirror- 
symmetric targets, and is not affected by external distur- 
bances inside the “clean” time window (i.e., before the 
arrival of reflections from the edge of the conducting 
plate). The shortcoming of the TDSRT is that it only 
applies to a smaller number of target types due to its 
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Figure 1. Measurement of Transient Electromagnetic 
Pulse Response Using TDSRT 


Key: 1. Monocone Antenna; 2. Monopole Antenna; 3. 
Conductive Plate; 4. Target; 5. Pulse Source; 6. Sampling 
Oscilloscope. 











limited measurement range (the linearity of a conducting 
plate is only a few meters). For example, it can measure 
a mirror-symmetric target, but not a tilted spheroid; it 
can measure vertical polarization (relative to the mirror 
image plane), but not horizontal polarization. 


By directly measuring target transient EM pulse response 
in free space, it is possible to remove the constraints 
imposed by field dimensions, target shape or polariza- 
tion. But to perform such measurement, a high- 
performance directional wideband pulse antenna is 
required, and many technical problems such as coupling 
between the wideband transmitting and receiving 
antennas, and elimination of the effects of background 
noise and interference must be solved. Therefore, in 
most experimental studies of target transient EM char- 
acteristics, the traditional TDSRT is generally used. 


Il. Measurement of Transient EM Pulse Response of 
the Target and the Antennas 


1. Transient Electromagnetic Measurement System 


The transient electromagnetic measurement system 
(TEMMS) we have developed’ uses the Tektronix TEK 
7904A/7S11/7T11 wideband sampling oscilloscope 
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(bandwidth: DC - 14 GHz) and PC-class microcom- 
puters. We have developed a waveform acquisition con- 
troller which can be directly plugged into the I/O slot of 
the microcomputer, a microcomputer-controlled target 
platform, a pulse generator, . wideband pulse antenna, 
and associated software packages. These components are 
integrated to form a multifunctional automatic measure- 
ment system (Figure 2) which is capable of measuring, 
collecting, processing, and performing spectral analysis 
of transient EM signals, as well as producing false-color 
images. This system has a vertical resolution of 12 bits 
and a selectable horizontal resolution of 8, 10, or 12 bits. 
We have constructed an experimental impulse radar 
based on the TEMMS and have used it to perform 
measurement and analysis of transient EM pulse 
response of a conducting sphere and a conducting pro- 
late spheroid in actual space (i.e., outer field); we have 
also derived the impulse response and calculated the 
wideband RCS of targets in the time domain. 


2. Design of Experiment 


The actual measurement environment is a complicated 
system which consists of the target, the background and 
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Figure 2. Block Diagram of the TEMMS 


Key: 1. Sampling unit; 2. Sampling and scanning unit; 3. Delay-line unit; 4. Sampling head; 5. Input; 6. Output; 7. 
External scan input; 8. Synchronous output; 9. Trigger input; 10. Receiving antenna; 11. Pulse generator, 12. 
Transmitting antenna; 13. Target; 14. Waveform acquisition controller; 15. Power supply for step-motor drive; 16. 
Step motor for target platform; 17. Hard disc drive card; 18. Monitor card; 19. Multifunctional card; 20. Hard disc; 
21. High-resolution color monitor, 22. Floppy drive unit; 23. Printer, 24. Plotter 
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Figure 3. Signal Flow Diagram of the Target, the Background and the Antenna 





the antenna. To extract the transient EM pulse response 
of the target from this system requires careful design of 
the experiment. 


(1) Measuring the Transient EM Pulse Response of the 
Target, the Background and the Antenna System 


The output signal from the pulse generator x, is radiated 
from the transmitting antenna and split into two parts: 
one part is radiated from the front of the antenna toward 
the target, and then scattered by the target into the 
receiving antenna; the other part is radiated from the side 
of the antenna and coupled into the receiving antenna 
(Figure 3a). The background signal x, enters the receiving 
antenna both from the front and from the side. Therefore, 
the received signal y (i.e., the EM response of the target, 
the background, and the antenna) also consists of three 
parts: the scattered signal from the target, the coupled 
signal from the side of the antenna, and the background 
signal; it can be expressed in the following form: 


y=X_*h,, *h,, *h*h,, *h, +X *h,, * hy *h, + X, 
* (h,, + h,,) (1) 


where * is the time-domain convolution operator, and 
all quantities are functions of the time t; h,, is the 
forward impulse response of the transmitting antenna, 
h,, is the side impulse response of the transmitting 
antenna, h,,, is the forward impulse response of the 
receiving antenna, h,, is the side impulse response of the 
receiving antenna, h,, is the side-path (space) impulse 
response between the transmitting antenna and the 
receiving antenna, h,, is the path (space) impulse 
response between the transmitting antenna and the 
target, h,, is the path (space) impulse response between 
the target and the receiving antenna, and h is the 
impulse response of the target. 


(2) Measuring the Transient EM Pulse Response of the 
Background and the Side Antenna System 


Part of the transmitting signal x, enters the side of the 
receiving antenna, and the background signal x, enters 
the receiving antenna from the front and from the side 
(Figure 3b). Therefore, the received signal (i.e., the 
transient EM pulse response of the background and the 


side antenna) y, consists of two parts: the antenna 
side-coupled signal and the background signal: 


Yo = Xo * hy * hy * h,, + x, * (h,, + h,,) (2) 


(3) Measuring the Transient EM Pulse Response of the 
Forward Antenna System 


The signal from the pulse generator x, is radiated from 
the transmitting antenna and enters the receiving 
antenna via the equivalent paths h,,, h,, (Figure 3c). The 
received signal (i.e., the transient EM pulse response of 
the forward antenna system) x can be expressed in the 
form: 


x=Xx)*h, *h,, *h, *h,, (3) 


tp 
(4) Signal Processing 


By subtracting equation (2) from equation (1) and 
applying the commutative law and associative law of 
time-domain convolution, one can obtain the transient 
EM pulse response of the target and the forward antenna 
system: 


Z=Y-Yo=X% *h,, *h,, *h*h,, * h,, = (Xp * hy, * hy, 
*h.*h,)*h 


rp 


Substituting equation (3) into the above equation gives: 


z= 7— nesen=[xt— t)h(sz)dr = Ah (4) 


— am {"e(t—1)( «dr 


is the time-domain convolution operator. The units of x 
and z are in mV, the unit of t is in ns, and the unit of h 
is in ns''. By measuring y, yo, and x respectively, one can 
obtain the impulse response of the target h by finding the 
inverse convolution from equation (4): 


h=[y-yo] * 'x=z*"'x (5) 


where * “' is the inverse convolution operator. 








3. Actual Measurement of the Transient EM Pulse 
Response 


The measurement system used in this experiment is the 
TEMMS. Both the transmitting and receiving antennas 
are continuous-loading, varying-slope horn antennas 
with horizontal polarization. The center of the transmit- 
ting antenna, the receiving antenna, or the target is 
approximately 5.3 m in height; the distance between the 
centers of two antennas is 3 m; the distance between the 
mid point of a line connecting the two antenna centers 
and the center of the target is 10 m. This setup provides 
approximate backscatter measurement. When the 
system is used to measure the transient EM pulse 
response of the forward antenna system, the distance 
between the transmitting and receiving antennas is 20.4 
m, which is approximately equal to the sum of the 
distances between the transmitting antenna and the 
target, and between the target and the receiving antenna, 
as shown in Figure 3a. A mercury gas discharge tube is 
used as the pulse generator.® Figure 4 shows the time- 
domain pulse waveform, which can be approximated by 
a Gaussian pulse. The test targets include a conducting 
sphere (diameter d = 60 cm) and a conducting prolate 
spheroid (axial ratio = 3, semi-major axis = 78 cm; the 
major axis is parallel to the ground, and the azimuth 
angle, i.e., the angle between the major axis and the 
incident ray, is chosen to be 90°, 45° and 0°, respec- 
tively). The target is placed on top of a non-metallic 
support rod of the target platform. The waveform acqui- 
sition process consists of taking 256 sample points of a 
waveform displayed on the sampling oscilloscope, and 
applying averaging to 100 consecutive waveform mea- 
surements in order to reduce the effects of system noise 
and background noise. Figure 5 shows the transient EM 
pulse response of the forward antenna system. Figure 9 
shows the transient EM pulse response of the two targets 
and the forward antenna system. 


III. Inversion of the Target Impulse Response 


1. Deconvolution Method 


Under the assumption of a time-invariant linear system, 
the target impulse response contains complete informa- 
tion about the target. From the impulse response one can 
calculate various target parameters such as RCS, reso- 
nance frequency and natural frequency. The target 
impulse response obtained from measured data plays an 
important role in target identification. 


One can see from equation (5) that the determination of 
target impulse response h from measured transient EM 
pulse response is a deconvolution problem. Because of 
noise in the measured data, the procedure of determining h 
by numerical deconvolution is unstable. There are several 
deconvolution methods: the Fourier transform method,’ 
the statistical characteristics method,® the conjugate gra- 
dient method (CGM),”'° the singular-value decomposition 
method,'' the time-domain filter method,’ the inverted 
spectrum deconvolution method’ and other methods." 
But each of these methods has certain shortcomings. 
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Figure 4. Pulse Source Waveform (attenuation 70 dB) 
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Figure 5. Transient EM Pulse Response of the Forward 
Antenna System 





In the CGM introduced by Tseng and Sarkar,'° the 
deconvolution problem is reduced to a problem of func- 
tion minimization. However, because of the difficulty in 
calculating the gradient of a high-dimensional vector 
function, it is approximated by only three terms. This 
so-called modified CGM is not a CGM in a rigorous 
sense. In 1985'° and 1987,° Tseng and Sarkar applied 
this method to determine the impulse response of a 
conducting sphere from measured transient EM pulse 
response data. The 1985 result was less than satisfactory; 
the basic characteristics showed an attenuated sine wave, 
which is quite different from the actual impulse response 
of a conducting sphere. The 1987 result showed strong 
oscillation in the tail of the waveform; this oscillation 
could be reduced by applying artificial weighting, but the 
result still differed considerably from the actual impulse 
response (see reference 13 or 14). 


2. Modified Steepest Descent Method (SDM) 


The modified SDM is an iterative method proposed by 
this author to determine the target impulse 
response;'*:'* its basic formulas are as follows: for an 
arbitrary time-domain function c = c(t) and corre- 
sponding vector c = [c(0), c(At), ..., c((N - 1)At)]’, At is 
the sampling interval of c(t) and N is the number of 
sample points. The definitions of the norm |] |] and A* 
are given below: 


el 


eof? = (c.c) = >) CeCnde) Fe 


eee [x6 — t)e(r)dt 
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The deconvolution problem of equation (5) is equivalent 
to the problem of finding the minir»um of the function 
F(h) = (Ah - z, Ah - z) =||Ah- Al’, where h = h(t), z = z(t), 
h = [h(0), h(At), ..., h((N - 1)At)]", z = [z(0), z(At), ..., z ((N 
- I)At}", Ah = [Ah),o, Abia, .-, Ablaviyat)” (samples 
have been taken of the time-continuous functions z and 
Ah). 


Having selected the initial value hy and specified the 
stopping rule, the following iteration procedure is car- 
ried out: 


2 -~ &—s.P.o RAR =) 14P [k= e + oP. 4 = 0.1.46) 


When the stopping rule is satisfied, then h,,, is the 
minimum solution desired (denoted by h), and the 
corresponding number of iterations is m = 1 + 1; other- 
wise let 1 — 1 + | and equation (6) 1s repeated. This 
function minimization approach is an approximate 
SDM in which the estimate of deconvolution error and 
the control of the number of iterations are based on the 
correlation coefficient.'*:'* For a finite-length discrete 
signal x(n), z(n) and h(n) (n = 0, 1, .... N - 1), numerical 
convolution can be performed: 


¥ 
: _ heey an Sela \Ele — 
2(a) = 2(n) h(a) = Dd z(arh(a—&), 


' 
- 


gs > O.lier — ] 


The iteration procedure is terminated when 


%-1 
r<esn . 


‘+ 
f i J 1 oo 
R, = p =i if > a(n ad Pp 
° L; $ ” l= rs | 


iL4 


where R... is a pre-selected threshold of the correlation 
coefficient. 


Numerical calculations show that selecting R,, = 0.98 gives 
satisfactory results, but better results are obtained by 
selecting R.. = 0.995. Increasing the value of R,, further 
will greatly increase the number of iterations. Therefore, in 
practice the choice R., = 0.995 is a reasonable compromise 
between achieving higher accuracy and limiting the 
number of iterations. If equation (7) is still not satisfied 
when m = 10, then the iteration procedure is terminated. 


3. Inversion of the Impulse Response of a Conducting 
Sphere and a Conducting Prolate Spheroid 


In Figure 6, the solid curve is the normalized impulse 
response of a conducting sphere obtained by applying 
inverse Fourier transform to the measured data in the 
frequency domain (the time origin is the time that corre- 
sponds to the maximum amplitude).'° The waveform of 
the impulse response of the conducting sphere consists of a 
negative impulse immediately followed by a positive step 
(contribution of the direct reflection), then followed by a 
negative pulse (contribution of the creeping wave). 


The measured waveform of the conducting sphere z 
(Figure 9a) and the antenna response x (Figure 5) are 
processed using the modified SDM (after 10 iterations) 
to obtain the impulse response h (Figure 10a). The 
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Figure 6. Impulse Response of a Conducting Sphere; 


solid line—h,,,,,,; dashed line—h,,.., 





correlation coefficient between z = x * h and z is R,. = 
0.98820. By comparing the normalized value of the 
inverted result, h,,.,..n (time is normalized with respect to 
t = d/c, where c 1s the speed of light in vacuum, and 
amplitude is normalized with respect to the maximum 
amplitude of the curve, as indicated by the dashed line in 
Figure 6) with h,,,,.,, one can see that the characte™istics 
of the two curves and the delay times of the creeping 
waves are essentially the same. 


Figures 10b-10d show the impulse responses of a con- 
ducting prolate spheroid (azimuth = 90°, 45°, and 0°) 
obtained by processing the measured waveform z (Fig- 
ures 9b-9d) and the antenna response x (Figure 5) using 
the modified SDM. The corresponding number of itera- 
tions m for each case is 8, 10, 10; the correlation 
coefficient R. for each case is 0.99636, 0.99382, 0.94753. 
The results corresponding to azimuth angles of 90° and 
45° are quite good, but the correlation coefficient at zero 
azimuth is rather low because in this case the nose of the 
spheroid is pointed at the antenna, hence the return 
signal is very weak compared to clutter and interference. 


It can be seen from Figures 10b-10d that the conducting 
prolate spheroid also exhibits the characteristics of direct 
reflection and creeping wave. As the prolate spheroid 
rotates from 90° to 0°, the direct reflection becomes 
weaker, hence the positive step gradually disappears, 
while the creeping wave changes from flat to sharp. The 
delay time of the creeping wave is determined by the 
path length of the creeping wave. The measured delay 
times corresponding to azimuth angles of 90°, 45° and 0° 
are respectively 8.3 ns, 10.2 ns and 10.8 ns, which 
compare favorably with the calculated delay times of 
7.54 ns, 9.69 ns and 11.01 ns. 


ee Measurement of Wideband Target 


1. Existing Problems With the Frequency-Domain 
Method of RCS Measurement 


There are two different methods of measuring target RCS: 
frequency-domain method and time-domain method. The 
traditional approach is to use frequency-domain method. 
In order to measure the wideband target RCS with this 
method, a sufficiently large number of frequency points 
are selected within the band, and single-frequency RCS 
measurements are obtained at successive points; these 











measurements are then combined to construct the wide- 
band target RCS. This method is similar to the frequency- 
domain method of measuring target impulse response; it is 
a complicated procedure that must be performed in a 
microwave anechoic chamber. 


2. Time-Domain Method of RCS Measurement 


If the frequency spectrum (i.e., the Fourier transform) of 
the target impulse h = h(t) is denoted by H = H(f), then 
the target RCS o can be expressed in terms of the 
amplitude spectrum | H] as follows:'® 


o = 4nr7]H|” (8) 


where r is the distance between the target and the 
receiving point. 


Let us define the normalized amplitude spectrum of the 
impulse response and the normalized RCS as follows: 


Horm = |H]/max{ H] (9) 
Snorm - _ (10) 


where the symbol max indicates the maximum value 
within a specified frequency range. 


By transmitting a narrow pulse toward the target and 
measuring the transient EM pulse response in the time 
domain, the wideband target RCS can be calculated 
using either the impulse response method or the tran- 
sient response method. 


(1) Impulse Response Method 


First, the target impulse response h is obtained by 
finding the deconvolution of the measured time-domain 
transient EM pulse response, as indicated in equation 
(5); then the amplitude spectrum |H] or the normalized 
amplitude spectrum H,,.,.,,, is determined; finally, the 
RCS o or the normalized RCS o,,.,,,,, is obtained from 
equation (8) or equation (10). But this method requires 
deconvolution of noisy data, which is an unstable pro- 
cess. For this reason, the impulse response method has 
not been widely used in practice. 


(2) Transient Response Method 


By applying a Fourier transform to equation (4), one 
obtains Z = X x H, where Z = Z(f), X = X(f) are 
respectively the Fourier transforms of the functions z = 
z(t) and x = x(t). It follows that 


H = Z/X, |H] =|ZjAX| (11) 


In the transient response method, first the transient 
EM pulse responses of the target and the forward 
antenna system z and of the antenna system x are 
measured; then z and x are Fourier transformed to 
obtain the amplitude spectra |Z] and |X], from which 
|H] can be determined using equation (11); finally, the 
RCS o is obtained from equation (8). For an actual 
time-domain waveform (particularly when noise is 
present), |Z] and |X] become very small as f increases, 
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therefore their relative errors become quite large, 
making the quotient of |Z] and |X] unstable. However, 
within the -20 dB band (f,-f,) of the amplitude spectra 
|Z] and |X] (i.e., when 


|Z|>0. 1 + max{/Z|},|X|>0.1 + max{/ X}} 


equation (11) can be used to provide a stable solution of 
the amplitude spectrum |H], from which the wideband 
RCS o can be determined. Thus, the transient response 
method is a preferred method for calculating wideband 
amplitude spectrum or RCS. 


Based on equation (9), the normalized amplitude spectra 
of Z and X are respectively 


Znorm = |Zj/max| Z] (12) 
X norm = | X|/max] X| (13) 


From equations (11)-(13), one obtains 
|H] ” B X ) eS (14) 


where 


B = max Z]/max| X] 


Using equation (14) and the definition of H.,,.,.,, from 
equation (9), we have: 


Horm - (Zncra! Nacemlnoren ( I 5) 


3. Determination of Target RCS Using the Transient 
Response Method 


Figure 7 shows the entire procedure of using the tran- 
sient response method to determine the normalized RCS 
of a conducting sphere from the measured transient EM 
pulse response. Figure 7a is the normalized amplitude 
spectrum of the transient EM pulse response z of the 
spherical target and the forward antenna system (Figure 
9a) (the spectrum is obtained using a 1024-point FFT; 
the -20 dB band is approximately 70-800 MHz; a small 
spectral peak at f = 0 is due to measurement error). 
Figure 7b is the normalized amplitude spectrum of the 
transient EM pulse response x of the forward antenna 
system (Figure 5) (the -20 dB band is approximately 
0-920 MHz). Figure 7c and Figure 7d are the normalized 
amplitude spectra of the conducting sphere obtained 
from equation (15), H,,,,, (at different selected fre- 
quency ranges). It can be seen from Figure 7c that above 
1 GHz large errors appear in the calculated spectrum due 
to the division operation between two small numbers. 
The resonance peak that appears in the vicinity of 1.14 
GHz apparently is a false signal caused by the fact that 
the normalized amplitude spectrum of the antenna 
system, X.orm (i.e., the denominator in equation (15), 
shown in Figure 7b), approaches zero at this frequency. 
But within the 70-800 MHz band, the calculated results 
of H,,.m are quite good (Figure 7d). Figure 7e shows the 
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Figure 7. Procedure for Determining the Normalized RCS of a Conducting Sphere Using the Transient Response 
Method 





normalized RCS o,,.,..,; the five resonance frequencies at 
the lower end of the spectrum are clearly visible in the 
figure. 


A comparison between the measured value and the 
theoretical value of the normalized RCS of a conducting 
sphere is shown in Figure 8. It can be seen from the 
figure that the transient response method produces 
smaller error in the normalized RCS than the impulse 
response method. For f below 0.5 GHz, the results of the 
transient response method are in complete agreement 
with theoretical results, whereas the results of the 
impulse response method clearly differ from the theoret- 
ical results. For the four resonance frequencies at the low 
end of the frequency spectrum (within the -20 dB band), 
the maximum relative errors in resonance frequencies 
between measured values and theoretical values for the 
impulse response method and the transient response 
method are 6.29 percent and 1.82 percent respectively. 
At higher frequencies, the measured values of o,,.,.,, tend 
to be too low because the conducting sphere used in this 
experiment is not an ideal sphere, and its relative losses 
are larger at higher frequencies. 


Figure 11 shows the normalized RCS of a conducting 
sphere and a conducting prolate spheroid obtained from 
measured data using the transient response method. It 
can be seen from Figures 11b-11d that as the azimuth 
angle of the prolate spheroid changes from 90° to 45° to 
0°, an increasing number of resonance frequencies begins 
to appear. In other words, the number of resonance 
frequencies increases with decreasing azimuth angle. 


V. Conclusion 


In this paper, a new technique for measuring and ana- 
lyzing the transient EM pulse response that is not subject 
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Figure 8. Comparison of the Measured and Theoretical 
Values of Normalized RCS of a Conducting Sphere 
Key: 1. Theoretical value; 2. Impulse response method; 

3. Transient response method 





to the limitations of the traditional TDSRT is presented. 
An experimental impulse radar which is based on the 
TEMMS is used to measure the transient EM pulse 
response of a conducting sphere and a conducting pro- 
late spheroid in actual space (i.e., outer field). 


A modified SDM is used to obtain the impulse responses 
of the conducting sphere and the prolate spheroid from 
measured transient EM pulse responses of the target and 
the forward antenna system and of the antenna system. 
The impulse response of the conducting sphere given in 
this paper is more accurate than previous results pub- 
lished in the literature; the impulse response of the 
conducting prolate spheroid is published for the first 
time. 
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Figure 10. Target Impulse Response 


Key: |. Conducting sphere; 
2. Conducting prolate spheroid 





Based on the measured time-domain waveforms of the 
transient EM pulse responses, the normalized amplitude 
spectra, the normalized RCS and the resonance frequen- 
cies of the conducting sphere and the conducting prolate 
spheroid are calculated using the impulse response 
method and the transient response method. The tran- 
sient response method introduced in this paper is an 
effective method for calculating wideband target RCS; it 
is also a practical engineering technique for calculating 
accurate target resonance frequencies. 
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[Text] 


Abstract 


The TEM horn for transmitting and receiving narrow 
pulse or impulse waveform has been discussed in this 
paper. A novel design concept and calculation method 
have been presented. According to this concept, a ridged 
conical TEM horn has been designed and fabricated. The 
experiment results show that its performance is excel- 
lent, so it is very valuable for the study of transient field 
and ultrawideband (UWB) radar. 


Key wor-ls: TEM horn, antenna, UWB radar, impulse 


1. Introduction 


TEM horn is an ultrawideband antenna which is widely 
applicated in the transmitting and receiving system of 
impulse waveforms. So far many theory analysis and 
design for this type of antenna have been presented,'~’ but 
they are limited to frequency domain. This is correct in 
some degree. For air medium, it is considered that TEM 
wave propagation is independent of frequency. By using 
Fourier transform, many characteristics such as transverse 
field distribution and wave impedence, can use straight- 
forwardly in time domain, but when discussing impedence 
match, it is inaccurate if we also use the concept of 
frequency domain. For example, many researchers think 
that the reflection is the least by using Chebyshev or 
Butterworth impedence transform in the internal of the 
horn, but they have ignored an important feature, namely 
time order of the narrow pulse. Above mentioned methods 
of impedence transform are constituted in the condition of 
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stationary sinusoidal! wave, all of the reflection wave at the 
impedence ladder join in the superposition procedure. 
However, the impulse wave will be reflected only when the 
pulse have arrived the impedence ladder. we should con- 
sider the time order of the superposition, generally 
speaking only a small section of the ladder takes part in 
superposition, so that above mentioned transform design 
has lost grounds in the condition of narrow pulse. Some 
authors think that the impedence of the horn aperture is 
the air wave impedence (1 20x), and take it as load impe- 
dence of the horn when they do impedence transform 
design. Sometimes this will bring on a disastrous result, a 
pulse train which 1s almost equal amplitude wil! be 
appeared. In fact, the match concept is not so simple, wave 
mode theory should be used for analysis strictly, as shown 
in Figure |. If we expand radiation wave of the horn in 
space spherical mode wave function,’ then we will notice 
that it is in terms of an infinite sets. and expressed as 


econ + P.E>) = 2b,E: 


2a,(H} — PLH>) = %,Hi _ 


in the above expression, the left indicate wave guide mode, 
superscript “+” means the external travelling wave, “-” 
means reflection wave, n is the mode ordinal number, a, is 
the nth mode amplitude, I, is its reflection coefficient. 
The right indicate the sum of space mode, and the symbol 
is as same as the left. The wave impedence for each mode 
not only relate to wave type, but also distance. In the very 
short distance, the impedence value of the low mode have 
closed to that of free space, the higher the wave mode is, 
the farther the required distance is. However, their wave 
impedence either tend towards zero or infinite on the short 
distance. Obviously, the wave impedence as a design value 
will relate to the size of the horn and its flared angle. 
Although it is difficult to compute it strictly, we should 
point out that the wave impedence of the horn will be more 
and more close to that of free space if the length of the horn 
is lengthened and its flared angle is enlarged, hence, the 
reflection of the aperture will be reduced. 


2. Theory 


The radiation waveform of a TEM horn ts required to 
have good fidelity, that is to say it is required to have 


OSD 
Q(Y0 


Figure 1. Internal and External Wave Mode Represen- 
tation of TEM Horn 
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little reflection. These reflections include that of the 
aperture, internal transmission, and feed terminal. It is 
difficult to calculate these reflections of the aperture, 
even though in frequency domain, there is not yet ready 
formula that can available. But we know that the 
reflection is little when the ratio of the size of the horn 
aperture to wavelength is large in frequency domain. It is 
not difficult to oo that there is similar case in time 
domain, when the horn aperture size is more large than 
the pulse width (to measure by light meter). TEM trans- 
mission line characteristic impedence is unmeaning for 
discussing the match between a horn and free space, for 
example, the characteristic impedence of two wire could 
be 377 Ohms when their distance is very close. Yet it is 
almost utter reflection. 


The reflection will be arised in the process of wave 
transmission in the horn. Here we should use the concept 
of characteristic impedence. For an impulse, the trans- 
mission line segment that wave arrived will have effect 
and the transmission line that the wave unarrived will be 
“invisible,” so that only a very short transmission line 
segment that wave have arrived will arise reflection for 
each instant, and the remain parts afterwards the seg- 
ment can be taken as a matching load, as shown in Figure 
2. If the reflection coefficient is constant, then the 
transmission coefficient is also constant, so the pulse 
waveform will be not distorted. This is the basic concept 
to discuss the impedence matching in time domain. 


We take Z,(0) which is the start terminal impedence of 
the TEM horn as the normalized value. Suppose the 
normalized characteristic impedence of transmission 
line is z(x), then the reflection coefficient of each point in 
the process of the wave march forward is 


_ 2(2 + Ar) — 2(2z) 
M2) = eadtdz)teG) 


Where z(x + Ax) is the average characteristic impedence 
of the segment that wave have arrived, 2Ax is the length 
of the segment. For the impulse waveform, the segment 
is more less than the all length of the horn, namely it is a 
minute-sized. So the change of the characteristic impe- 
dence of the segment can be regarded as linear, and their 
average value is namely the middle value. From it we 


have 
tA 
1 \dz 
M2) = > s(x) 





that is 
z(2x) on gis . (3) 


where a = 2/Ax. 


Thus it can be seen that if let the I'(x) equals a constant, 
then the TEM horn should be exponentially flared, and 
when I(x) = 0 the TEM horn characteristic impedence is 
constant. This is more simple than the wideband impe- 
dence matching in frequency domain. 


a(x) |ie(z+Ar) 
1242! 





Figure 2. Discontinuity Equivalent Circuit of the TEM 
Hom 





There are many methods, such as variation method and 
conformal mapping method can be used to compute the 
characteristic impedence of the TEM horn. For a com- 
plex shape horn, it is the most convenient and reliable to 
use the integral equation method. Figure 3 shows the 
cross section of a ridged conical TEM horn, its top plate 
and low plate are equal potential surface, the potential 
equation can be written as 


9,8) = [ G@.? y0@ yas (4) 


where 7’ is source point, 7 is field point, p is the 
density of line charge, G.(%, 7) is Geree 
function of potential 


PY —F 
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ee 








(5) 


where « is the permittivity, 7” is the mirror point 
of 7 against the x axis, 











Figure 3. Cross Section Coordinate of Ridged Conical 
TEM Horn 





In the expression (4), p is unknown function. Since the 
conductor is equal potential, we can let the 


Ar,0)=1, ons, 
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thus integral equation can be written as 


[ce ~ 7) pds = 1, ons (6) 


The moment method can be used to solve integral 
equation into a matrix equation, and the charge density 
p(r’[—-+ over r]) can be evaluated, so the total charge will 
be obtained by integration on s 


Q= | ow )ds’ (7) 


The distribution capacity of the horn is 
c=Q/2 (8) 
Thus characteristic impedence can be found as 
z= I/ev = 2/Qv (9) 
where v is the velocity of light. 


Figure 4 is the charge distribution on the antenna 
computed by moment method, the radius of the horn is 
normalized as 100. A quarter of the horn has been shown 
in the figure. 


The abscissa represents the direction of x axis, ab is 
thickness of the ridge, bc is highness of the ridge, d the 
fringe point. Figure 5 shows the characteristic impe- 
dence of the horn under the case of different highness of 
the ridge, the abscissa is highness of the ridge. 


Since only a few of points can be computed by numerical 
method, but the realistic structure should be continual 
geometry curve, for it we can build interpolation formula 
by using Lagrange method or Newton method to imple- 
ment its technical design. 


3. TEM Horn With Large Flare Angle 


We have designed and fabricated an ultrawideband 
TEM horn with large flare angle for use as a feed element 
in impulse radar system. It is made by posted a copper 
foil on the glass fibre reinforced plastic. A photograph is 
shown in Figure 6. The diameter of the horn aperture is 
50 centimeters, and its length is 43 centimeters, the flare 
angle is 60 degrees. A transformer from coaxial to double 
line has been used to feed the horn. A ridge which is 
exponentially flared has been adopted to guide success- 
fully the TEM wave into the horn. 
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Figure 4. Charge Distribution of the TEM 
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Figure 5. Characteristic Impedence Under Different 
Ridge Width 
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Figure 6. Photograph of the Ridged Conical TEM Horn 





A large flared angle and large aperture horn transmits 
TEM wave to have obvious special feature, and it 
appeared the geometry optical characteristic in fre- 
quency domain, namely the radiation field is mainly 
determined by the field of the aperture corresponding 
point on the range of the horn flared angle, this point is 
exactly that the ray pass the point. This property still be 
kept in time domain, namely the impulse waveform on a 
point in space is determined by the waveform of corre- 
sponding point of the horn aperture. Because the wave- 
form of the horn aperture is basically undistortion, thus 
the waveform for all direction of the horn flared angle is 
basically unchanged. This is its distinguishing feature. 


The experiment results show that the impedence of the 
horn has good performance. With the exception of a few 
frequencies, the input VSWR of the antenna is 2:1 from 
0.3 GHz to 7 GHz, its Smith chart is shown in Figure 7. 
Their time domain waveforms shown as in Figure 8. In 
this experiment, a 5 centimeter high conical antenna 
used as the transmitting antenna, the TEM horn is used 
as receiving antenna. The source pulse is Gauss wave- 
form with 150 ps rise time and 350 ps pulsewidth (-3 
dB), bottom width is 500 ps. Since effect of differential, 
receiving waveforms have positive and negative peaks, 
total width approximate | ns. It is seen from the wave- 
forms in time domain that the waveform shape and 
width is basically undistortion and its trail is very small, 
hence the fidelity of the waveform is excellent. 


4. Conclusions 


TEM horn is an ultrawideband antenna. So far almost all 
of the study for it is used the concept of frequency 
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Figure 7. VSWR Plot for 0.3 GHz to 7 GHz for the 
TEM Horn 





domain, and only limited to the plate TEM horn. So a 
novel design concept and computation method, which is 
based on the wave shape fidelity, is presented in this 
paper. According to the theory, charge distribution and 
characteristic impedence for the TEM horn have been 
calculated by numerical method, and designed and fab- 
ricated a new type of conical ridge TEM horn antenna. 
The experiment results show that this horn has good 
impedence performance and excellent waveform 
fidelity. This is very valuable for the study of transient 
field and ultrawideband radar. 
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[Abstract] Subaperture range-Doppler (R-D) processing, 
an improvement over conventional R-D processing, is 
an effective fast (6 minutes on AST-286 microcomputer) 
high-accuracy algorithm in microwave imaging, one of 
the premier techniques for radar target signature mea- 
surement. The basic principles are introduced; then 
computer simulations are done and the data compared 
with experimental results. High, well-distributed resolu- 
tion and uniformity of amplitude have been obtained. In 
the computer simulation, reflectivity amplitudes of -15 
dBsm are obtained. The experimental testing was done 
at the BUAA Electromagnetic Engineering Laboratory, 
equipped with a model BH-205I RCS/2-D imaging 
system; the test object was a YB-47B scale model (1/28 
size). 


Figures 1-2, not reproduced, are coordinate systems for 
parametric measurement. Figures 3-5 are reproduced 
below, as is Table 1, listing imaging data. 
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Figure 3. Target Imaging Results for Six Ideal Points 
via Subaperture R-D Imaging. 
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Figure 4. Taget Imaging Results for YB-47 Scale Model 
via Subaperture R-D Imaging. 
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Figure 5. Target Imaging Results for YB-47 Scale 
Model via Conventional R-D Imaging. 











x(cm) Y(cm) RCS (dBsm) 
. pl 0. 78 0. 78 —15. 19 
p2 3. 94 3. 94 — 15. 12 
p3 40. 16 0. 78 —15. 00 
p4 35. 4 3. 94 — 15. 08 
ps —7. 87 40. 16 —15. 03 
pé6 3. 94 35. 4 —15. 01 
Table 1. Microwave Imaging Data; x=horizontal coordi- 
nate, y=vertical coordinate. 
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Stealth 


Enhanced Aircraft Angular Glint: New Aircraft 
Stealth Technique 


94P60269A Beijing YVHANG XUEBAO [JOURNAL 
OF ASTRONAUTICS] in Chinese Vol 15 No 2, Apr 94 
pp 80-87 


[Article by Yin Hongcheng [3009 4767 2052] of the 
University of Electronic Science & Technology of China 
(UESTC) and Huang Peikang [7806 1014 1660] of the 
Second Academy, China Aerospace Industry Corp. 
(CASC): “Enhanced Aircraft Angular Glint Method: A 
New Aircraft Stealth Technique”; MS received 22 Oct 92] 


[Abstract] A general 2-D problem of electromagnetic 
scattering from a multiple cavity-backed longitudinally 
loaded slotted perfectly conducting cylinder with arbi- 
trary cross section is formulated. The formulation is used 
to account for the scattered field from scale models, for 
which a series of experiments has been conducted in an 
anechoic chamber. Excellent agreement between theory 
and experiment is obtained. Analysis of the experimental 
results near resonant frequency shows that the reduction 
of backscattering radar cross section (RCS) over wide 
aspect angles and obvious enhancement of target angular 
glint (linear deviation of angular position in radar target 
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line of sight) can be achieved by impedance loading such 
as multiple cavity-backed slots to realize a new aircraft 
stealth technique. Experimental model fabrication and 
microwave testing were aided by colleagues at the Bei- 
jing Institute of Environmental Features. 


Figures 1-2, not reproduced, show field computations by 
a moments method. Figure 3-7 are reproduced below. 
The authors conclude that from measured results for a 
TE-mode 12.35-GHz incident wave and a U-type cavity 
slotted loaded wall cylinder scale model (equivalent to 
aircraft wing loading), within a +/- 20 deg attitude-angle 
range, the RCS reduction exceeds 20 dB, with a max- 
imum point reduction of 35 dB. Also for this same 
incident-angle range, the probability that normalized 
angular-glint linear deviation e bar exceeds | (full target 
size) is 59.2 percent, and maximum e bar is 15.2, a value 
permitting realization of an effective new stealth tech- 
nique. There are no tables. 
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3. Harrington, R.F., ‘““Time-Harmonic Electromagnetic 
Fields,” New York: McGraw-Hill, 1961. 














Figure 3. Scale Model of U-Type Cavity Slotted Loaded Wall Cylinder (Aircraft Wing); cylinder length=200 
(all units are in mm). 
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Figure 4. Comparison of Backscattering Theoretically Calculated (dotted line) and Experimentally Measured (solid 
line) for f=12.75 GHz TE-mode Incident Wave on Double-Slotted Cylinder and Solid Cylinder. 


Key: (a). RCS, with o‘=incident angle, o(,) = RCS in dBsm; (b). Normalized linear deviation e bar. 
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Figure 5. Comparison of Backscattering Theoretically Calculated (dotted line) and Experimentally Measured (solid 
line) for f=9 GHz TE-mode Incident Wave on U-Type Cavity Slotted-Wall Cylinder and Solid Cylinder. 


Key: (a). RCS; (b). Normalized linear deviation. 
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Figure 6. TE-Wave (f=8.75 GHz) Backscattering Experimental Data for Double-Slotted Cylinder and Solid Cylinder. 


Key: (a). RCS; (b). Normalized linear deviation e bar. 
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Figure 7. TE-Wave (f=12.35 GHz) Backscattering Experimental Data for U-Type Cavity Slotted-Wall Cylinder and 
Solid-Wall Cylinder. 


Key: (a). RCS; (b). Normalized linear deviation e bar. 
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New Method for Electromagnetic Performance 
Design of Structural Absorbers 


94P60273A Chengdu DIANZI KEJI DAXUE XUEBAO 
[JOURNAL OF UNIVERSITY OF ELECTRONIC 
SCIENCE AND TECHNOLOGY OF CHINA 
(UESTC)] in Chinese Vol 23 No 1, Feb 94 pp 53-56 


[Article by Gao Chunyang [7559 2504 7122], Zhao Bolin 
[6392 0130 3829}, and Rao Kejin [7437 0344 6210} of 
the Institute of Information Materials Engineering, 
UESTC, Chengdu 610054: “New Method for Electrical 
Performance Design of Structural Absorbent Materials”; 
MS received 7 Jul 93, revised 18 Aug 93] 
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[Abstract] A new method for designing EM performance 
of structural absorbers is presented. Unlike a conven- 
tional CAD process, the new method is based on a 
Microsoft Fortran analysis program which simulta- 
neously calculates reflection and transmission coeffi- 
cients and reasonably arranges impedance transforma- 
tions between layers of a multilayer structure to ensure 
wideband (2-20 GHz) operation. Design examples using 
existing materials—domestically made hybrid absorb- 
ers—exhibit performance much better than usually 
obtained via a CAD process. With the new method, 
reflectivity R and transmissivity T of 3-or-4-layer struc- 
tural materials are both under -10 dB, and unidirectional 
material thickness is under 20mm. The structural 
absorbers include glass-fiber-reinforced plastic or Kevlar 
fiber (first layer). 


Figure | below shows a 3-layer structural absorbent 
material with standard EM parameters and thicknesses 
labeled, while Figure 2 below shows a 17-mm-thick 
4-layer design and Figure 3 shows a 3-mm-thick 3-layer 
(not including the thick back lining, or fourth layer) 
thin-wall absorber used for wing front and back edges, 
air intakes, and other lip openings. 














Figure 1 3-Layer Structural Absorbing Material 





N=4 


i. =1.03-j0.02 
p,.=1.27-j 0.65 
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ra =1.5-j 08 
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Z .,=0.521 7 exp( —0.566 3) 
Z ,,=0.505 7 exp (0.9020) 
Z ,; =0.387 7 exp ( -—0.170 8) 


Z ,«=0.2757 exp ( -0.7537) 


Figure 2. A 4-Layer Design. 
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ec oo uw, =1.03-j0.02 ¢,,=4.02-j04 d,;=1 mm 

-2 8 9 10 ti 12 13 14 15 16 17 18 19 w,,=15-jOS ¢,,=45-j13 d,=1 mm 

J/ GHz u.s=1$-j072 €,,=130-j62 d,;=1 mm 

Figure 3. A 3-Layer Thin-Wall Design. 

Reference: and a number of double- layer dielectric specimens have 


1. Zhao Bolin, Gao Zhengping, Rao Kejin, “Reflection 
and Transmission Coefficients for Electromagnetic 
Waves Obliquely Incident on Layered Lossy Mediu- 
m”[{in Chinese], DIANZI KEJI DAXUE XUEBAO 
[JOURNAL OF UESTC], 1993, 22(5):483-489. 


Bandpass Radome Research Highlighted 


Measurement of Complex Permittivity of 
Double-Layer Dielectrics 
94P60281A Beijing TONGXIN XUEBAO [JOURNAL 
OF CHINA INSTITUTE OF COMMUNICATIONS] 
in Chinese Vol 15 No I, Jan 94 pp 46-58 


[Article by Xia Jun [1115 6511] and Liang Changhong 
[2733 2490 3163] of Xidian University, Xian 710071: 
“Measurement of Complex Permittivity of Double- 
Layer Dielectrics by Electromagnetic Open Resonator,” 
supported by grant from State Education Commission’s 
Ph.D. Fund; MS received 21 Jun 91, revised Apr 92] 


[Abstract] A new technique for measuring the complex 
permittivity of double-layer dielectric specimens—such 
as are used on radomes, radar absorbing materials 
(RAM) for other aircraft surfaces, and microwave inte- 
grated-circuit substrates—at microwave and millimeter- 
wave frequencies using an electromagnetic (EM) open 
resonator is proposed. At the 8-mm band, an EM open 
resonator automatic measurement system has been built, 


been measured. Finally, the error is theoretically ana- 
lyzed and some calculated examples are given. 


Figure | below shows the structure of the loaded open 
civity, while Figure 2 below shows the intra-cavity 
dielectric interface approximation curved surface. 
Figure 9 below shows a schematic of the 8-mm measure- 
ment system. Figures 3-8 and 10-16, not reproduced, 
show the following: a perturbation volume element, an 
equivalent small volume, two graphs of dielectric- 
specimen thickness modification vs t (thickness of outer 
dielectric layer), two graphs of dielectric-specimen thick- 
ness modification vs other EM parameters, four more 
graphs of various EM parameters, and three schematics 
of the actual measurement method, respectively. Table | 
below lists measured results. 
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Figure 1. Schematic of Loaded Open Cavity. 
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Figure 9. Schematic of 8mm-Band EM Open Resonator M 


Key: 1. Printer; 2. Monitor, 3. HP8341B frequency synthesizer, 4. HP8349B amplifier, 5. HP83554A frequency- 
doubling module; 6. QBJF-01 isolator, 7. WL501 level-meter; 8. SL7151 digital voltmeter, 9. MT60 grating length 
gauge; 10. IAAS digital display; 11. BCD-GPIB [binary coded decimal-general purpose inter-face bus] converter. 
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113. 9077 57766 4.3867 | 0.000041 





115. 2322 37836 4. 3662 0. 000069 





116. 4050 21000 4. 6207 0. 000136 





116. 0045 28060 4. 4631 0. 000052 





115. 4977 14193 4. 4389 C. 000078 





























116.8832 | 104245 4.4682 | 0.000060 


Table 1. Measured Results. 


Key: 1. High-pressure polyethylene; 2. Polytetrafluoroethylene; 3. microcrystal glass. 


Resonant frequency=35.0000 


GHz, Y-axis mode number q=27, cavity resonator length = 117.1817mm, cavity quality factor Q,=120323. 
*System-measured results for single-layer speciment. ** Cavity resonator length and quality factor of cavity after 


loading with double-layer dielctric specimens. 
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Simulation of Antenna Pattern by Complex Rays 
94P60281B Beijing TONGXIN XUEBAO [JOURNAL 
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[Article by Ruan Yingzheng [7086 4481 6927] of the 
Dept. of Microwave Engineering, UESTC, Chengdu 
610054: “Simulation of Antenna Pattern by Complex 
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91] 


[Abstract] Based on complex ray theory, the Gaussian 
beams emitted by complex source points are used to 
simulate the far-field paiicrn of antennas. Three 
methods of simulation are proposed—i.e. mainlobe 
matching, peak matching, and point-by-point match- 
ing—and some calculated results for several varieties of 
patterns are given. Numerical results of analysis and 
optimization show that the peak matching method is an 
accurate and practical approach for antenna pattern 
simulation. 


Three figures, not reproduced, show a complex- 
source-point simulated pattern via mainlobe matching 
(with aperture length L = 5 x wavelength), the same via 
peak matching, and the same via point-by- point 
matching, respectively. Table 1, not reproduced, lists 
typical antenna patterns and the error introduced by 
complex-ray approximation of these patterns. 
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